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ABSTRACT 

We present the full catalog of Young Stellar Objects (YSOs) identified in the 18 molecular clouds 
surveyed by the Spitzer Space Telescope “cores to disks” (c2d) and “Gould Belt” (GB) Legacy surveys. 

Using standard techniques developed by the c2d project, we identify 3239 candidate YSOs in the 18 
clouds, 2966 of which survive visual inspection and form our final catalog of YSOs in the Gould Belt. 

We compile extinction corrected SEDs for all 2966 YSOs and calculate and tabulate the infrared 
spectral index, bolometric luminosity, and bolometric temperature for each object. We find that 326 
(11%), 210 (7%), 1248 (42%), and 1182 (40%) are classified as Class O+I, Flat-spectrum, Class II, and 
Class III, respectively, and show that the Class III sample suAFers from an overall contamination rate by 
background AGB stars between 25% and 90%. Adopting standard assumptions, we derive durations 
of 0.40 — 0.78 Myr for Class 0-hI YSOs and 0.26 — 0.50 Myr for Flat-spectrum YSOs, where the ranges 
encompass uncertainties in the adopted assumptions. Including information from (sub)millimeter 
wavelengths, one-third of the Class 0-hI sample is classified as Class 0, leading to durations of 0.13—0.26 
Myr (Class 0) and 0.27 — 0.52 Myr (Class I). We revisit infrared color-color diagrams used in the 
literature to classify YSOs and propose minor revisions to classification boundaries in these diagrams. 

Finally, we show that the bolometric temperature is a poor discriminator between Class II and Class 
III YSOs. 

Subject headings: infrared: stars ISM: clouds stars: formation - stars: low-mass 


1. INTRODUCTION 

Most of the current star formation within 500 pc of 
the Sun occurs in the Gould Belt, a ring of nearby O- 
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type stars inclined approxiniately 20° with resp ect to 
the Galactic Plane (jHCTSche]l Il847t iGoiildl Il879f l . The 


Sun is located interior to the Gould Belt, which has a ra- 


dius of about 350 pc le.g..|Glubelll967HStothers & Froeel 

1974 

JComeron et al.lll992t Ide Zeeuw et al.lll999HP6DDe] 

2001 

). All of the nearby, well-studied giant molecu- 


lar clouds (e.g., Ophiuchus, Taurus, Perseus, Serpens, 
Orion, etc.) are located with in this ring (see Figure 1 of 
IWard-Thoiupson et al.l[2?)??l for a schematic diagram of 
the Gould Belt in galactic coordinates). 

Due to the high visual extinction toward giant molec¬ 
ular clouds, a complete census of star formation in these 
clouds requires infrared observations. The Infrared As¬ 
tronomical Satellite (IRAS) provided the first such un¬ 
biased survey of the Gould Belt (and the rest of the 
sky), but was limited in both sensitivity and spatial res¬ 
olution. For young stellar objects (YSOs) still in the 
embedded phases of evolution, IRAS was only sensitive 
to objects with lu minosities above 0.1 (d/140pc)^ Lq 
( jMvers et al.lll98'ii1 l. For more evolved YSOs, IRAS was 
even less sensitive. The spatial resolution provided by 
IRAS ranged from 0.5' at 12 ^m (0.02 — 0.07 pc at dis¬ 
tances of 150 — 500 AU) to 2' at 100 ^m (0.09 — 0.29 pc at 
distances of 150 — 500 AU), insufficient to resolve individ¬ 
ual stars forming in clustered regions of molecular clouds. 
Ground-based near-infrared surv eys (e.g., the Two Mi- 
cron All Sky Survey (2MASS); iSkrutskie et al.l I200 (tI 1 
provide much higher spatial resolution but are unable 
to detect the most deeply embedded, heavily extincted, 
and/or lowest luminosity YSOs. Ground-based mid- 
infrared studies are generally limited by the atmosphere 
to only the brightest and most luminous objects. Finally, 
selected regions of the Gould Belt were mapped in the 
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mid-infrared by the Infrared Space Observatory (ISO). 
While ISO surveys identified many new deeply embed¬ 
ded and low luminosity YSOs, the y only focused on small 
portions of molecular clouds (e.g. jBontemps et alJ 120011 : 
IPersi et ^1120031 : IKaas et aI1l200^ . 

The launc h of the Spitzer Space Telescope 
(|Werner et alJ l2004ll in 2003 opened new avenues 
for infrared surveys of giant molecular clouds. With 
its combination of high sensitivity and large format 
detector arrays operating at 3.6 — 160 /rm, Spitzer was 
the first mid-infrared facility capable of mapping entire 
molecular clouds within reasonable amounts of time. 
The Spitzer Legacy project “From Molecular Cores to 
Planet-forming Disks” (Cores to Disks, or c2d) took 
advantage of these features by mapping seven large, 
nearby Gould Belt molecular clouds and nearly 100 
iso lated dense cores. An overview of the survey is given 
bv lEvans et al.l (1200311 and a summary of res ults for the 
large clouds is given by lEvans et al.l (1200911 . The c2d 
survey identified over 1000 YSOs, many of which were 
newly identified, and demonstrated Spitzer’s ability to 
detect YSOs that are too deeply embedded, too low 
luminosity, and/or too low mass to be detected by 
previous facilities. 

The Spitzer Gould Belt (GB) Legacy survey followed 
up the c2d survey by mapping an additional 11 Gould 
Belt clouds using the same observation and data reduc¬ 
tion strategies as c2d. Combined, the two surveys tar¬ 
geted nearly every major site of star formation within 500 
pc of the Sun. The only two exceptions are the Taurus 
and Orion Molecular Clouds, both of which are so large 
they were the s ubject of their own , dedicated surveys b y 
different teams (|Rebull et al.ll201(il : lMegeath et al.ll2012ll . 
Adding the results of the GB survey to those obtained 
by c2d nearly triples the number of YSOs identified in 
the Gould Belt. A preliminary attempt to combine the 
results from the c2d, GB, Taurus, an d Orion surveys 
was presented bv iDunham et all (l2014f) . but a full com¬ 
bination, including reconciliation of different methods for 
observing, data reduction, source extraction and classifi¬ 
cation, is left for future work that combines the expertise 
of the respective surveys. Here we only consider the c2d 
and GB surveys. 

While the full list of YSOs identified by Spitzer in 
the c2d and GB surveys has been used in several 
studies, including th ose on clustering of young stars 
(iBressert et al.l 1201011, star f ormation rates and thresh¬ 
olds ( Heiderm an et a l1l201C)f) . the luminosities of proto¬ 
stars JDunha^^^lT]2^^j_j;ehitionships between gas 
and star formation (|Evans et al.ll20l4) . and the physi- 
cal nature of empirically d efined classes of young stars 
((Heiderman fc EvansI 120151) , this list itself has not yet 
appeared in publication. To rectify this, we present in 
this paper the full list of YSOs in the combined c2d and 
GB clouds. For the c2d clouds, this list updates and re - 
places the previous version given bv lEvans et al.l (1200911 . 
including an analysis of contamination by background 
AGB stars, updated extinction corrections, and revised 
spectral energy distributions (SEDs). Although we use 
this YSO list to investigate classification of young stars 
and the durations of these classes, the primary purpose 
of this work is to present a full catalog of Spitzer c2d+GB 
YSOs. While useful on its own, this YSO list will reach 
its full potential in combination with results from the 


far-infrared and submillimeter surveys of the Gould Belt 
by the Herschel Space Observatory and SCUBA-2 on 
the James Cle rk Maxwell Telescope, once they become 
fully available (fWard-Thomoson et al.ll2007tlAndre et al.l 
120101: iHarvev et al.ll2013l) . Bv providing the full Spitzer 
YSO list, we enable such comparisons. 

The organization of this paper is as follows. We de¬ 
scribe our method in (J2l including an overview of the 
c2d and GB surveys (' (12.11) . the procedure for identifying 
YSOs (' (12.2|) . residual contamin ation in our YSO sam¬ 
ple by background AGB stars 1 (12.31) . and the construc¬ 
tion of full SEDs and corrections for extinction 1 (12.41) . 
We present our results in (J31 including the full list of 
YSOs 1(13.11) and a discussion of classification of young 
stars 1(13:^1331). We discuss timescales for the classes 
of young stars in (jd) Finally, we summarize our results 
in ^ 


2. METHOD 

A concise summary of these surveys, the adopted pro¬ 
cedure for identifying young stellar objects, and the 
methods used to construct SEDs and to correct the pho- 
to metry for foreground extinction are given in Section 2 
of IDunham et al.l (120131 ). Here we repeat the key infor¬ 
mation and refer the reader to Dunham et al. for more 
details. 


2.1. Overview of the Surveys 

The Spitzer c2d survey (PI: N. J. Evans) imaged 
seven large, nearby molecular clouds in the Gould Belt 
and approximately 100 isolated dense molecular cores. 
lEvans et al.l (I2003D review both the primary science 
objectives o f c2d and its basic observation plan, and 
lEvans et al.l (|2009D summarize the results in the large 
clouds. The Spitzer GB survey (PI: L. E. Allen) is a 
follow-up program that imaged an additional 11 molec¬ 
ular clouds in the Gould Belt, completing most of the 
remaining clouds in the Gould Belt except for the Tau¬ 
rus and Orion molecular clouds (see ©• Both sur¬ 
veys imaged their respective targets at 3.6- 8.0 pm with 
the S pitzer Infrared Array Gamera (IRAG; lEazio et al.l 
120041) ■ and at 24-160 yum i mages with the Multiband 
Imaging Photometer fMlPS: lRieke et aT1l2004D . A data 
pipeline developed by c2d was applied to the observations 
obtained by both programs. This pipeline includes data 
reduction, source extraction, and band-merging the re¬ 
sults from each wavelength to produce final images and 
source catalogs, and has been described in detail else¬ 
where (jHarvey et al.ll200^ lEvans et al.ll2007D . 

For each cloud targeted by c2d or GB, Table [1] lists 
the parent survey, the assumed distance and reference for 
the distance, the total area mapped, the total cloud gas 
mass, and references to in-depth studies of each cloud 
based on the c2d or GB data, when available. The 
Spitzer Astronomical Observation Requests (AORs) for 
each cloud can be found in these references to individual 
cloud studies; Appendix A lists the AORs for the clouds 
without published data references. The isolated cores 
surveyed by c2d are not considered here except for a few 
cases where the cores are actually parts of larger cloud 
complexes and included in the results for those clouds 
(namely Cepheus an d Ophiuchus North: iKirk et al1l2009l : 
iHatchell et al.l[2012D . Note that the distances to these 
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clouds are not all well determined, and some cloud dis¬ 
tances are still under debate. One such example is the 
uncertainty over the distance to Aquila, with two dis¬ 
tances (260 pc an d 429 pc) commonly adopted in the lit¬ 
erature (see( e.g.. iHarvev et al.ll2 0n6| iGutermuth et al.l 
120081: iDzib et ^120111 : iMaurv et al.M201l[ for detailsl. 
For this reason we do not list formal distance uncertain¬ 
ties and instead refer the reader to the references given 
in Table [T] for full discussions of the distances to each 
cloud. 

The total cloud areas mapped and total cloud gas 
masses contained with in these areas are taken from 
iHeiderman et al.l (I2010D . with the masses determined 
from extinction maps (see Heiderman et al. for details). 
For most clouds, the mapped areas were chosen to map 
the entire cloud above a minimum visual extinction in 
the range 2 < Ay < 3, with the exact value for each 
cloud depending on the total area that could be mapped 
within the allocated time. For Serpens and Aquila, only 
the cloud areas above Ay = 6 were mapped due to time 
constraints and difficulty in separating the cloud from 
the general extinction close to the Galactic midplane. 


2.2. Identification of Young Stellar Objects 

The data reduction pipeline used by both the c2d and 
GB surveys creates band-merged source catalogs where 
the sources extracted at each Spitzer wavelength are 
matched with each other and with the 2MASS cata¬ 
log, creating a final catalog of 1.25 — 70 /rm photome¬ 
try for each cloud. Can didate YSOs are ident ified using 
a method developed bv iHarvev et al.l (l2007all and sum¬ 
marized in all of the publications focused on individual 
clouds that are listed in Table [T] In this method, the 
Spitzer SWIR E Legacy survey of t he ELAIS N1 extra- 
galactic field (|Lonsdale et al.l[2003D is processed to sim¬ 
ulate the same sensitivity and extinction distribution of 
the c2d and GB clouds, and then used to locate back¬ 
ground galaxies observed through the target clouds in 
various infrared color-magnitude diagrams. Using these 
results to define the expected colors and magnitudes of 
background galaxies, each source extracted in the cloud 
maps with infrared colors that can not be explained as 
an extincted background star is assigned an unnormal¬ 
ized “probability” of being a galaxy or YSO based on its 
position in each color-magnitude diagram, its morphol¬ 
ogy in the two shortest Spitzer IRAC bands (3.6 and 
4.5 /im), and its 24 and 70 /im flux densities. Since it 
is unnormalized, this quantity is not a true probability 
but a relative parameter expressing the likelihood of be¬ 
ing a background galaxy (high er values indicate hig her 
likelihood of being galaxies; see IHarvev et al.l[2007ar for 
details). The final boundary between candidate YSO 
and candidate galaxy in this quantity is set to provide a 
nearly complete elimination of SWIRE sources, likely at 
the cost of exclu ding YSOs with ov erlapping colors and 
magnitudes le.g.. iHsieh fc Laill2013[l . 

One potential drawback to the c2d-|-GB method of 
identifying YSOs is that it requires detections in all four 
IRAC bands (3.6 — 8.0 ^m) as well as the first MIPS 
band (24 /im). As the 24 /im MIPS band has lower res¬ 
olution and sensitivity, some of the fai ntest and/or most 
clust ered YSOs may be missed (e.g., IGutermuth et al.l 
120091 Gutermuth et ah, in preparation.). We note that 
several alternative classification methods for separating 


background galaxie s and YSOs have been presented in 
the literature (e.g.. IGutermuth et al.ll2009HRebull et al.l 
120101 : iKrvukova et al.1 120121 : iHsieh fc Tail 1201311 . While 
they are all based on a similar philosophy of separating 
sources in various colors and magnitudes, they differ in 
their detailed implementations of this philosophy. A full 
comparison of these methods and ultimate consensus on 
best practices for extracting YSOs from infrared surveys 
of star-forming clouds is left for future work. 

Applying these selection criteria results in the identi¬ 
fication of 3239 candidate YSOs in the 18 c2d and GB 
clouds. We note here that we only consider sources that 
are located in the overlap regions imaged by both the 
IRAC and MIPS instruments, as these are the sources 
that can be most reliably classified using the method 
summarized above. For technical reasons, the MIPS ob¬ 
servations often covered much larger areas (see the ref¬ 
erences listed in Table [1] for more details on the obser¬ 
vations of each individual cloud). Several of the indi¬ 
vidual cloud studies listed in Table [T] also present candi¬ 
date YSO lists in the MIPS-only regions using alterna¬ 
tive identification and classification techniques, but such 
sources are not considered here. We visually inspected 
each of the 3239 candidate YSOs to remove resolved 
galaxies a nd image artifacts that are sometimes misiden- 
tified fsee lEvans et al.l[2009l for details). Additionally, a 
few known YSOs that were saturated or undetected at 
one or more wavelengths were missing from the initial 
list and were added by hand. In the end, we were left 
with a final list of 2966 YS Os, and since all ha ve passed 
visual inspection, we follow lEvans et all (I2009D and drop 
the prefix “candidate” at this point. 

Finally, we note that initial results from Herschel sur¬ 
veys of the c2d and GB clouds suggest that a mod¬ 
est number of the most deeply embedded protostars are 
undetected by Spitzer, increasing t he number of Class 
0 protostars^° by 15% - 30% (e.g.. IHarvev et al.l [20131 : 
iStutz et aT]l2013t iSadavov et al.ll2014ll . Revisions to the 
YSO lists presented here will thus be necessary once com¬ 
plete Herschel results are available for each cloud. 

2.3. Background Giant Star Contamination 

Our final sample of YSOs is likely contaminated 
by background gia nt stars with infrared excesses. 
lOliveira et al.l (|2009ll performed optical spectroscopy to¬ 
ward 78 YSOs in Serpens and identified 20 (26%) as back¬ 
ground Asymptotic Giant Branch (AGB) stars. While 
the overall contamination rate is relatively low. Figure [I] 
shows that the contamination fraction strongly depends 
on the infrared spectral index Starting from pos¬ 

itive values, as a decreases, the contamination fraction 
remains close to zero until the Class H/HI boundary is 
reached at a = —1.6 (see ^3.2l for a full definition and dis¬ 
cussion of YSO classes), and then steeply rises to >50% 
for more negative values of a. On average, 62% of the 
Class HI sources examined by Oliveira et al. are actually 

Class 0 protostars are the youngest class of YSOs; see ^S.ll for 
the formal definition of this class of objects. 

The quantity a is defined as the slope of the infrared SED 
in log(AS'A) vs. lo g(A) and is used to classify YSOs, as discussed 
in detail in and Ea While a calculated from extinction 

corrected photometry is used in later sections, here we use the 
values calculated from the observed photometry for consistency 
with the previous studies to which we compare. 
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Figure 1 . Left: Histograms of the infrared spectral index (a, uncorrected for extinction) for the YSOs observed bv lOliveira et all 1I20Q9I) . 
The shaded histogram shows the full sample of 78 objects whereas the dashed histogram shows the 20 objects identified as background 
AGB stars by Oliveira et al. Right: The black symbols and line plot the contamination fraction in each a bin, where this fra ction is defined 
as the number of background AGB stars in this bin divided by the total number of objects in this bin, using the results from lOliveira et alJ 
l|2QQ9lj . The colored symbols plot the contamination fraction for each Glass (red for Glass O+I, green for flat-spectrum, blue for Class II, 
and purple for Class III). The Class II and Flat-spectrum results are plotted at the mean value of a for each class, whereas the Class III 
and Class O-l-I results are plotted at the mean values of the boundaries for each Class and the display ranges of the x axis. The dotted line 
in each panel marks the Class III boundary (q < —1.6). 


background AGB stars, whereas only 5% of the Class II 
sources they examined were found to be contaminants. 

Very high rates of AGB star contamination amon g 
Class III sources are also found bv lRomero et al.l (j20I2l l. 
also with optical spectra. They found Class III contami¬ 
nation rates of 100% based on observations of 30 objects 
in Lupus V and VI. While the fact that they only ob¬ 
served objects that met initial selection criteria for tran¬ 
sition disks may bias their results, they observed 45% 
and 39% of the total Lupus V and VI Class III popula¬ 
tions, respectively; thus the true contamination fractions 
are at least ^40% in each cloud. 

The true contamination fraction is likely a strong 
function of galactic latitude, and the specific clouds 
considered above (Serpens, Lupus V, and Lupus VI) 
all lie within 2 — 1 0° of the galactic plane. Indeed, 
iRomero et al.l (j2012t l found a lower contamination rate 
of 25% for Ophiuchus {b = 18°), consistent with the 
Ophiuchus contam ination rate of 20% measured by 
iCieza et'all (|2010ll using similar methods. To further 
investigate the contamination rate in our full sample of 
YSOs, Figure [2] plots a Spitzer color-color diagram con¬ 
sisting of [3.6]-[24] versus [3.6]-[4.5], using the observed 
photometry. Inspection of Figure [5] shows that the Class 
III YSO population can be divided into a population 
of objects clustered at lower (bluer) values of both col¬ 
ors, and a second, more distributed population extending 
to higher (redder) color s. Using op t ical sp e ctra of can¬ 
didate tran sition disks. ICieza et (j20ini l. ICieza et'H] 
(|2012ll . and iRomero et al.l (1201211 f^nd that background 
AGB stars are predominately found at low values of [3.6]- 


[24], suggesting that the clustered population in Figure 
[2] may be dominated by AGB stars whereas the more 
distributed population may be dominated by YSOs. In¬ 
deed, the above studies found that all of their objects 
with [3.6]-[24] < 1.5 were in fact AGB stars, with de¬ 
creasing contamination rates at higher colors. We find 
that 298 of the objects in our final YSO sample have 
[3.6]-[24] < 1.5, all but three of which are classified as 
Class III objects. This represents 25% of the Class III 
YSO sample and 10% of the total YSO sample. We con¬ 
sider these objects as likely AGB stars and mark them as 
such throughout the remainder of this paper. However, 
we do not remove them from the sample since the studies 
cited above only targeted a subset of objects that were 
identified as transition disk candidates and thus there is 
no guarantee that their finding of 100% contamination at 
the lowest values of [3.6]-[24] applies to our full sample. 
Confirmation of the status of these objects requires ei¬ 
ther spectroscopic or proper motion studies. Such work 
should be considered as a high priority for future investi¬ 
gations seeking to refine the YSO sample presented here. 

The 25% Class III contamination by AGB stars found 
above is almost certainly a lower limit to the full number, 
since some conta minating AGB stars are found to have 
[3.6]- [24] > 1.5 (|Cieza et al.l l20ld . I2012t iRomero et al.l 
1201211 . However, these studies target samples that are too 
small and too biased to determine contamination rates as 
a function of [3.6]-[24] color. To obtain an upper limit to 
the possible contamination fraction. Figure [3] compares 
the distributions of a for our YSOs with those identi¬ 
fied in young stellar clusters bv iGutermuth et all (1200911 . 
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Figure 2. [3.6]-[24] versus [3.6]-[4.5] for all 2966 objects in our 

final YSO sample, using the observed photometry. The color of 
each symbol shows its spectral class deter mine d usi ng t he extinc¬ 
tion corrected infrared spectral index (see ani and E2] for details 
on classification), with red indicating Class 0-hI, green indicating 
Flat-spectrum, blue indicating Class II, and purple indicating Class 
III. The horizontal dashed line indicates the boundary of [3.6]- 
[24] = 1.5 used to identify highly likely AGB stars contaminating 
our YSO sample. The vertical line of sources at [3.6]-[4.5] 0.5 is 

a rounding artifact that is discussed in detail in Appendix [BI 

using a selection method that is also based on infrared 
colors and magnitudes but different in its details (see 
Gutermuth et al. for details). The c2d+GB sample con¬ 
tains a much higher fraction of Class III objects. While a 
full comparison of different methods for identifying YSOs 
is beyond the scope of this work, under the assumptions 
that the Gutermuth et al. method perfectly separates 
out background AGB stars and our regions contain the 
same fractions of sources in each Class, 87% of the Class 
III objects in the c2d -|-GB sample wou l d need to be re¬ 
moved to match the iGutermuth et al.l (|2009ll Class III 
fraction. We note here that the Gutermuth et al. sur¬ 
vey targeted smaller regions centered on dense clusters 
whereas the c2d and GB surveys mapped much larger 
areas, including those more distant from cluster centers. 
As the Class II and III YSO populations are generally 
more distributed than earlier Classes, the intrinsic frac¬ 
tion of Class III sou rces is likely high e r in th e c2d and 
GB surveys than in IGutermuth et al.l (1200911 . For this 
reason, our above assumption that the surveys contain 
the same fractions of sources in each Class results in an 
upper limit to the Class III contamination fraction. 

Synthesizing all of the above information, our sample 
of Class III YSOs suffers from between 25% and 90% con¬ 
tamination by AGB stars. These values are lower and up¬ 
per limits, respectively, for the reasons listed above, and 
the most likely value is somewhere in between. Tighter 
constraints on the true contamination rate will require 
large optical spectroscopic campaigns and/or proper mo¬ 
tion studies coupled with revisions to the selection cri¬ 
teria used to identify YSOs in infrared surveys. These 


Figure 3. Histogram showing the distribution of infrared spectral 
index for the YSOs identi fied in the c2d-|-GB clou ds (this work, 
shaded histogram) and by IGutermuth et a D 120091) in a survey of 
young stellar clusters (dashed histogram). The dotted line marks 
the Class III boundary (o < —1.6). 

contamination rates must be kept in mind when consider¬ 
ing the results presented here and in other studies based 
on the c2d-|-GB YSO list (such as those discussed in (JT]). 

2.4. Constructing SEDs and Correcting for Extinction 

Foll owing the method first summarized bvlEvans et al.) 
(|2009ll and updated by iDnnham et al.l (|2ni,‘I[l . we com¬ 
pile as complete an SED as possible for each of the 2966 
YSOs. We refer to Dunham et al. for a full description 
of this process, but give a summary here. To construct 
the SEDs, we included: 

1. Selected ground-based optical and infrared pho¬ 
tometry as compiled by the authors of the indi¬ 
vidual cloud papers listed in Table [T] 

2. Wide-field Infrare d Survey Explorer (WISE; 
iWright et al.l 1201(1( 1 12 and 22 /rm photometry 
from the all-sky catalog^^. 

3. Spitzer 160 /rm photometry for sources that are 
detected but not located in saturated or confused 
regions. Since the c2d pipeline does not extract 160 
/rm sources, we measured flux densities using aper¬ 
ture photometry (using the IDL procedure APEr) 
and standard aperture corrections as given by the 
MIPS Instrument Handbook^^. 

4. Selected 350 /rm photometry from a tar geted sur¬ 
vey of protostars with SHARC-lP'^ (|Wu et all 
I2007L A. Suresh et al. 2015, submitted). 


Available at http://irsa.ipac.caltech.edu/cgi-bin/Gator/nph-scan7mission 
Available at http://irsa.ipac.caltech.edu/data/SPITZER/docs/ 
mips / mipsinstrumenthandbook / 

The Submillimeter High Angular Resolution Camera II 
(SHARC-II) was a 350 /rm b olometer array oper ated at the Caltech 
Submillimeter Observatory IIDowell et al.|[2b03l) . 
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5. Selected 450 and 850 /im photometry from 
the S CUBA^^ Legacy Catalog (|Di Francesco et ^ 

[200l . 

6. Other submillimeter and millimeter photometry 
from the literature, where availabl e from large 
surveys of nearby molecular clouds (lYoung et ah 
[20M[20MlEnoch et al.ll2006l [200^ iBelloche et al. 
l2011bllal:lXurTetal.ll2nilh 


Regarding the last three items, we have access to 
complete submillimeter or millimeter surveys for only 
6 out of the 18 clouds (Chamaeleon I, Chamaeleon 
II, Chamaeleon III, Ophiuchus, Perseus, and Serpens), 
plus a partial survey of Aquila and piecemeal cover- 
age of other clouds from the SCUBA Legacy Catalog 
dDL Francesc o et al.||2 008[ l and our own 350 /rm observa¬ 
tions ( Wu et al.ll2007l Suresh et al. 2015, submitted). 
We consider a Spitzer source to match a submillimeter 
or millimeter detection if it is located within one beam of 
the detection, using the appropriate beam size for each 
survey. If multiple Spitzer sources match a single sub¬ 
millimeter or millimeter detection, we do not attempt to 
split the flux density and simply match it to the closest 
Spitzer source unless the SED clearly indicates a better 
match is obtained with a different Spitzer source. 

Finally, we correct all photometry for fo reground ex¬ 
tinctio n following the procedure outlined by lEvans et al.l 
(|2009f) . We only wish to correct for the foregound extinc¬ 
tion from both the clouds and the material between us 
and the clouds, and not the local extinction from dense 
cores surrounding protostars (this extincted emission is 
reprocessed to longer wavelengths and thus included in 
our compiled SEDs). Thus, we adopt extinction values 
as follows: 


1. Whenever possible, we adopt extinction values 
from the literature for Class II and III YSOs (classi¬ 
fied via infrared spectral index; see ^3.21) . as deter¬ 
mined from optical or near-infrared spectroscopic 
studies of each cloud that provide reliable spectral 
type (and thus extinction) measurements. We re¬ 
fer to the data references listed in Table [T] for more 
details on these studies. 

2. For all Class II and III YSOs with no published ex¬ 
tinction values, we de-redden to the intrinsic near- 
infrared colors of an assumed spectral type of K7. 
This particular spectral type is chosen because it is 
found t o be representative of the YSOs in the c2d 
clouds (jOliveira et al.l[200M 120101 see also Evans et 
al. 2009 for details). 

3. We calculate the mean extinction toward all Class 
II YSOs in each cloud, and then adopt that mean 
value for each of the Class O-l-I and Flat-spectrum 
YSOs in that cloud. 


With extinctions assigned as a bove, we de-redden the 
SED of each YSO using the iWeingartner fc Drain^ 
(|2001ll extinction law for Ry = 5.5. Whereas 

The Submillimeter Common-User Bolometer Array (SCUBA) 
was a 450 and 850 /im bolometer array operated at the James Clerk 
Maxwell Telescope. 


Evans et'aLl (l2009fl adopted a version of the 
Weingartner fc Draind (|200lD Ry = 5.5 extinction 
law that only included dust absorption, here we update 
the procedure to correct for both dust absorption and 
scattering. The choice of this particular extinction law 
is motivated by results showing that it is appropriate 
for the dense molecul ar clouds in which stars form (e.g., 
iChaoman et al.l[2009D . While we do caution that our 
approach only provides approximate extinction correc¬ 
tions, especially for the Class O-l-I and Flat-spectrum 
sources where the corrections only account for the mean 
cloud extinction, it is nevertheless the best that can 
currently be done and is significantly more reliable than 
ignoring the effects of extinction altogether. 

3. RESULTS 

3.1. Full List of Young Stellar Objects 

For each of the 2966 YSOs, Table H] lists a running in¬ 
dex, the cloud in which the YSO is found, the Spitzer 
source name (which also gives the source coordinates), 
the visual extinction used for extinction corrections, and 
the infrared spectral index (a), the bolometric tempera¬ 
ture (Tboi), and the bolometric luminosity (Lho/)- These 
last three quantities are calculated using both the ob¬ 
served and extinction corrected photometry, and are de¬ 
noted with primes in the latter case (a', Tboi', and Lboi')- 
The second-to-last column of Table [2] indicates whether 
or not each object has [3.6]-[24] < 1.5 and is thus likely 
a background AGB star, as described above in >12.31 

The final column of Table [H indicates whether or not 
each object is associated with a dense core as traced 
by submillimeter or millimeter continuum emission at 
A > 350 pm. As a lack of such an association does 
not necessarily imply that no core is present (see be¬ 
low), to avoid any confusion in the interpretation of this 
column we list either Y or no entry rather than Y or N. 
iDunham et al.l (|2013ll used such associations to define 
the sample of protostars in the c2d-|-GB clouds, with 
the rationale behind this definition being that detections 
in existing submillimeter and millimeter surveys of star¬ 
forming regions trace dense cores but not circumstellar 
disks due to the relatively low spatial resolutions and 
mass sensitivities of these surveys (see iDunham et al.l 
120131 for details). By requiring a dense core but making 
no assumptions about the infrared colors of protostars, 
this definition ensures a reliable sample of protostars that 
includes those viewed edge-on down outflow cavities, but 
comes at the cost of incompleteness to those protostars 
that are embedded in low-mass cores or located in clouds 
with incomplete or no avail able submillimeter or m illime- 
ter continuum surveys Isee IDunham et ani20l4 for de¬ 
tails). Associations with dense cores and the resulting 
sample of protostars must be revisited once the JGMT 
SGUBA-2 survey of the Gould Belt is complete, the first 
resul ts of which are ju st starting to become available 
(e.g., iPattle et al.ll201^ . 

The infrared spectral index (a, a') is defined as the 
slope of the infrared SED in log(AS'A) vs. log(A), where 
A is the wavelength and S\ is the flux density at that 
wavelength: 


d\og{XSx) 

dlog(A) 


( 1 ) 
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Extinction Corrected Bolometric Temperature (K) 


Figure 4. The extinction corrected bolometric luminosity plotted versus the extinction corrected bolometric temperature for 

all 2966 YSOs. The color of each symbol shows its spectral class determined from a' (the extinction corrected infrared spectral index), with 
red indicating Class 0+1, green indicating Flat-spectrum, blue indicating Class II, and p urple indicating Class III. Black symbols indicate 
sources with [3.6]-[24] < 1.5 that are considered likely to be background AGB stars (see ^2.3|l . Filled circles mark sources associated with 
submillimeter or millimeter continuum emission tracing d ense cores (see t e xt for details), while the plus signs mark sources with no such 
associations. The three thick lines plot model tracks from lYoung &: Evand (I20Q5I) for singular isothermal spheres of different initial masses 
undergoing inside-out collapse with 100% efficiency (all mass initially in the cores end up in the stars). These models end when infall stops 
and the cores hav e fully accret e d ont o the stars and are thus only showing evolution in the protostellar stage. The three thin lines plot 
model tracks from iMvers et ^ d 19981) for cores of different initial masses collapsing with accretion rates that decre ase exponentially w ith 
time and efficiencies less than one (i.e., only a fraction of the mass initially in the cores ends up in the stars, see IMvers et al.lll998l . for 
details). As above, these mod els are only relevant for the pr otostellar stage of evolution. The heavy da shed line on the le ft is the ZAMS 
from 0.1 to 2 Mq taken from i^Antona &: Mazzitellil 1119941) . The class boundaries as defined by IChen et al.l dl995l) are plotted as 
vertical dashed lines. 


In practice, we calculate a with a linear least squares fit 
to all available 2MASS and Spitzer photometry between 
2 and 24 /rm. The bolometric temperature of a source 
is defined to be the temperature of a blackbody with 
the s ame flux-weighted mean frequency (| Myers fc LaddI 
|1993[) . It is in essence a protostellar equivalent of stellar 
effective temperature; it starts at very low values 10 
K) for cold, starless cores and eventually increases to the 
effective temperatnre of the central (proto)star once the 
core and disk have fully dissipated. Lf^oh Li,oi , and 
Ti,oi are calculated by integrating over the full SEDs: 


Lboi = 47rd^ 



Thai = 1.25 X 10-“ 


/p°° vSvdv 

r s.d. 


K 


( 2 ) 

(3) 


In practice, we calculate these integrals using the trape¬ 
zoid rule to integrate over the finitely sampled SEDs. 

Figure ID plots Ll „, versus (a “BLT” diagram; 

I Mvers fc LaddI [TOOM for all 2966 YSOs identified in the 
c2d and Gould Belt clouds. The luminosities of Class 
0 and Class I sources are distributed over greater than 
4 orders of magnitude and extend to much lower lumi- 
nosites than predicted by simple evolutionary models. A 
full discussion o f this “protostellar lumi n osity problem” 
can be foun d inlOffner fc McKe ^ (1201111 . iDuiiham et al.l 
(IMl . and iDunham et al.l (|2014fl . but to summarize, 
mass-dependent and/or time variable mass accretion is 
required to explain the observed luminosities of proto¬ 
stars. Class II and Class III objects exhibit both a very 
narrow range of bolometric temperatures and a bifurca¬ 
tion in T/gj into two groups. Both effects are artifacts in¬ 
troduced by our extinction corrections and are discussed 









































in Appendix [Cj To summarize, the first effect is because 
the bolometric temperature of a Class II or III source 
depends more on the spectral type of the object than 
the amount and characteristics of circumstellar material, 
and for most objects a constant spectral type of K7 is 
assumed (see above). The second effect (bifurcation) is 
a consequence of optical data only being available for a 
subset of our YSOs. 

Tables [3] and |4] list, for each YSO, the same running 
index as in Table O followed by flux density and flux 
density uncertainty pairs for the 2MASS (1.25, 1.65, and 
2.17 /im) and Spitzer (3.6, 4.5, 5.8, 8.0, 24, and 70 p,m) 
wavelengths. All flux densities and flux density uncer¬ 
tainties are listed in mJy and rounded to two significant 
digits. Table |3] lists the observed values while Table |4] 
lists the extinction corrected values. The observed and 
extinction corrected flux densities that we have compiled 
at all other wavelengths in order to construct SEDS, as 
described above, are tabulated in Appendix ID] 

3.2. Classification of Young Stellar Objects 

Over the past few decades a standard model for low- 
ma ss star format i on has emerged. This model, developed 
bv lAdams et al.l (|1987[ 1 and summarized by iShu et aH 
(j 19871) . resulted from the merger of an empirical clas¬ 
sification system based on the infrared spectral indices 
of YSOs (iLada fc Wilkind[l984ll with theoretical st udies 
of the stages of the collaps e of a dense, rotating core (IShul 
Il977t l^rebev et ^119841 1. Recent reviews on classifica¬ 
tion of young stars and the associations between observed 
SE P classes and phy sical evolutionary stages ar e given 
by lEvans et al.l (l200^ and iDunham et al.l (1201411 

Including the revisions bv iGreene et al.l (1199411 . these 
empirical Classes are as follows: 

• Class 0-fI: a' > 0.3; 

• Elat-spectrum: —0.3 < a' < 0.3; 

• Class II: —1.6 < a' < —0.3; 

• Class III: a' < —1.6. 

We denote the infrared spectral index as a' to emphasize 
that we only consider extinction corrected values in this 
study, and we refer to the first SEP class as Class O-l-I 
instead of the more commonly used Class I to empha¬ 
size that Class 0 and I objects are both included in this 
definition (see below). 

Table [S] lists, for each cloud, the numbers of YSOs in 
each Class using the extinction corrected infrared spec¬ 
tral indices, both with and without the likely background 
ACB stars removed. Summed over all clouds, out of the 
2966 total YSOs, 326 (11%) are classified as Class 0-fI, 
210 (7%) are classified as Flat-spectrum, 1248 (42%) are 
classified as Class II, and 1182 (40%) are classified as 
Class III. If we first remove the 298 sources identified as 
likely background ACB stars, out of the 2668 remaining 
YSOs, 326 (12%) are classified as Class O-l-I, 210 (8%) 
are classified as Flat-spectrum, 1245 (47%) are classified 
as Class II, and 887 (33%) are classified as Class III. We 
emphasize that there are large variations between clouds 
in the relative numbers in each Class. To reinforce this 
point visually. Figures [5] and [6] plot, for each cloud, four- 
color bars where each color corresponds to a Class and 


All YSOs 


Auriga/CMC 
Chamaeleon III 
IC5146 


Ophiuchus North 




Perseus 




Cepheus 



Corona Australis 




Serpens 




Ophiuchus 



Lupus IV 


Chamaeleon II 
Lupus I 
Mu sea 
Aquila 

Chamaeleon I 
Lupus III 
Lupus VI 
Lupus V 
Total 


Figure 5. A four-color bar is plotted for each cloud, where each 
color corresponds to an evolutionary Class (similar to Figure [JJ 
red corresponds to Class O+I, green corresponds to Flat-spectrum, 
blue corresponds to Class II, and purple corresponds to Class III). 
The fractional length of each color represents the fraction of YSOs 
in each cloud with the corresponding Class. The clouds are sorted 
by Class O-hl fraction, with the largest fractions at the top and the 
smallest fractions at the bottom, and all 2966 YSOs are included. 

the fractional length of each color represents the fraction 
of YSOs with the corresponding Class. Figure [S] displays 
the results including all 2966 YSOs whereas Figure [6] 
displays the results after removing the 298 likely back¬ 
ground ACB stars. Finally, we also emphasize that our 
method of YSO identification is only capable of identi¬ 
fying Class III sources with detectable infrared excesses 
and is thus incomplete to the full populations of pre-main 
sequence stars in these clouds. 

As noted above, we refer to the first Class as Class 
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Excluding Likely AGB Stars 


Chamaeleon III 
Auriga/CMC 
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Figure 6. Same as Figure^ except after remo ving t he 298 objects 
identified as likely background AGB stars (see ^2.3|l . 

O+I, but th e original definition of this Class was sim - 
ply Class I (|Lada &: Wilking|ll984 iGreene et al.lll99ll . 
Class 0 objects were later added as an earlier Class for 
protostars too deeply embedded to detect in the in¬ 
frared bu t inferred through other means (i.e., outflow 
presence; lAndre et al.l Il993ll . They are defined as ob¬ 
jects with a ratio of submillimeter to b olometric lu¬ 
mino sity {Lsmm/Lboi) greater than 0.5% (lAndre et al.l 
[HH, where Lsmm is calculated for A > 350 fj,m. While 
they are defined observationally, the intended corre¬ 
sponding physical definition of Class 0 objects are pro¬ 
tostars with at least 50% of their total mass still re¬ 
sidin g in the surrounding, infalling core (|Andre et al.l 
1199311 . An alternative observational definition is given 
in terms of the bolometric temperature as objects with 



Extinction Corrected Bolometric Temperature (K) 

Figure 7. Extinction corrected infrared spectral index (a') plot¬ 
ted versus extinction corrected bolometric temperature for 

all 2966 YSOs, with the colors and symbols the same as in Fig¬ 
ure |4l T he horizontal l ines s how the class boundaries in a' as 
revised by IGreene et al.1 ll 19941) and adopted in this work, and the 
vertical dashed lin es show the Ti,oi class boundaries as defined by 
IChen et al.l II1995I 1. 

Tbol < 70 K (|Mvers fc LaddI [19931 : IChen et all 1199511 . 
While Lsmm/Lboi is generally a better discriminator be¬ 
tween Class 0 and Class I objects than Tbni, and is less 
sensitive to source geomet ry (e.g.. lYoung fc Evansll2005l : 
iDunham et al.llToifll . 120141 , it cannot be calculated with¬ 
out 350 /im photometry. As we do not have 350 /rm pho¬ 
tometry available for many of our YSOs, we only consider 
Tbol in this study. 

With the sensitivity of Spitzer, Class 0 sources are eas¬ 
ily d ete cted in the i nf rared fe. g.. INor iega-Cr espo et all 
2004 lYoung etlT 12004 IJorgensen et alJ 120061 : 
Dunham et al. " EoOt , 120081 : iTobin et al.l 1201011 . Fig¬ 
ure [7] plots a' versus for all 2966 YSOs. In general, 
the two quantities are anti-correlated, although there 
are some notable exceptions. There are a small number 
of SED Class O-bl sources (plotted in red in Figure [7|) 
with sufficiently high values to be classified as Class 
11 or 111 by None of these objects are associated 

with dense cores and are likely simply more evolved 
YSOs viewed through large extinctions that are not fully 
extinction-corrected (recall that, for SED Class 0-1-1 
sources, only the average cloud extinction is accounted 
for). There are also a few SED Class 11 YSOs (plotted 
in blue in Figure El with sufficiently low values 
to be classified as Class 0 by T^^oi- These are likely 
true protostars viewed at nearl y pole-on inclinations 
through outflow cavi ti es Ce.g.. IWhitnev et al.l 120031 : 
IRobitaille et al.l 120061 : ICrapsi et al.l 1200*^ ^though 
detailed follow-up study is required to rule out chance 
alignments between more evolved YSOs and dense cores. 

While there is a general anti-correlation between a' 
and Tjg;, this trend weakens for objects classified as Class 
O-bl by a' (objects with a' > 0.3). With substantial over- 
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lap in a' for objects classified as Class 0 or Class I by 
it is clear that these two Classes overlap in their infrared 
spectral indi ces. These results confi rm those presented 
previously bv lEnoch et al.l (|2n09[l and lEvans et al.l (j20n9[ l 
and indicate that mid-infrared data alone are insufficient 
to separate Class 0 and Class I objects; instead, full SEDs 
are required. 

Using our compiled SEDs, and restricting ourselves 
only to those YSOs associated with a dense core (both 
to restrict ourselves to confirmed protostars and to en¬ 
sure available submillimeter or millimeter detections for 
reliable measurements of our measurements 

result in 65 Class 0 protostars and 129 Class I pro¬ 
tostars. Thus, we find that approximately one third 
(65/[65 -I- 129] = 0.335) of protostars are classified as 
Class 0, in agreement wit h previous results based on the 
small er c2d-only sample (lEnoch et al.ll200^ lEvans et al.l 
I2009f) . While future work must revisit this ratio once 
complete, sensitive submillimeter coverage is available 
for all clouds through the Herschel and JCMT Gould 
Belt surveys, preliminary Herschel results suggest only 
modest increases t o the number of Class 0 objects (gen¬ 
erally 15% - 30%. iHarvev et al.l[20T^ IStutz et all 120131 
ISadavov et al.ll2014ll . leading to only small upward revi¬ 

sions in the ratio of Class 0 to Class I protostars. 

3.3. Infrared Color-Color Diagrams 

Figure |8] presents an extinction corrected infrared 
color-color diagram of all YSOs identified in this 
work constructed from the four IRAC wavelengths: 
[3.6]' — [4.5]' plotted versus [5.8]' — [8.0]'. The differ¬ 
ent Classes of YSOs (classified via a') are ge nerally well- 
separa ted in this diagram, as first noted by lAllen et al.l 
(|2004ll . In cases where photometry is only available at 
3.6 — 8.0 /im and thus the infrared spectral index can¬ 
not be reliably calculated, this diagr am can be used to 
approximately classify YSOs. Indeed, lAllen et all (120041 1 
used the following two conditions to define a box where 
Class II sources are located: 

1. 0.4 < [5.8] - [8.0] < I.l, and 

2. 0.0 < [3.6] - [4.5] < 0.8. 

The origin a l def inition of this “Class II box” by 
lAllen et ^ (120041) used observed photometry. Inspec¬ 
tion of Figure |5] reveals that minor revisions are required 
when using extinction corrected photometry. While we 
acknowledge that polygons with greater than four sides 
could possibly give better results, for simplicity we only 
consider rectangular regions. To define a new Class II 
box, we allow all four boundaries to vary and select those 
which simultaneously maximize the probability that a 
Class II YSO is located within the box (defined as the 
number of Class II YSOs located within the box divided 
by the total number of Class II YSOs) and the proba¬ 
bility that a YSO located within the box is classified as 
Class II (defined as the number of Class II YSOs located 
within the box divided by the total number of YSOs lo¬ 
cated within the box). In practice, we implement these 
definitions by selecting the boundaries which maximize 
the sum of these two probabilities. Our proposed new 
Class II box, plotted in panel (b) of Figure El has the 
following boundaries: 


1. 0.45 < [5.8]' - [8.0]' < 1.55, and 

2. 0.00 < [3.6]' - [4.5]' < 0.70. 

We note that the shape of th is revised bo x is qu alitatively 
similar to that proposed bv iGutermuthI (|2005D . who fol¬ 
lowed a similar procedure. With these revised bound¬ 
aries, the fraction of Class II YSOs located within the 
box increases from 84.5% to 89.9%, and the fraction of 
YSOs in the box that are Class II increases from 78.8% 
to 82.4%. 

Figure[9]presents a second extinction corrected infrared 
color-color diagram of all YSOs which uses both the 
IRAC and MIPS observations by plotting [3.6]' — [5.8]' 
versus [8.0]' — [24]'. Again the different Classes of YSOs 
(classified via a') are generally well-separated, as first 
noted bv iMuzerolle et all (I2004D . Indeed, Muzerolle et 
al. defined three regions in this color-color diagram to 
separate Class Ill/stellar, Class II, and Class O-l-I YSOs; 
these regions are marked in panel (a) of Figure El and 
are defined in Table m As above, Muzerolle et al. used 
observed photometry, and we find that minor revisions 
are required when using extinction corrected photometry. 
Furthermore, the original boundaries defined by Muze¬ 
rolle et al. did not include a region for Flat-spectrum 
objects. Thus, we follow the same procedure as above, 
again only considering rectangular regions, and maximize 
the sum of the eight probabilities (two for each Class), 
subject to the following constraints: (I) the minimum x 
and y Class III boundaries be held fixed at the Muzerolle 
et al. values, (2) there are no maximum x and y Class 
O-l-I boundaries, and (3) the regions can not overlap. Our 
revised boundaries are plotted in panel (b) of Figure E] 
and listed in Table El Table El also lists each of the two 
probabilities considered for both the original and revised 
boundaries. 

Ultimately, no region in either color-color diagram will 
classify YSOs with a perfect one-to-one correspondence 
to the infrared spectral index, but in cases where in¬ 
sufficient photometry at 2 and/or 24 /rm is available, 
approximate YSO classification is possible via infrared 
color-color diagrams. 

4. TIMESCALES FOR YOUNG STELLAR OBJECTS 
4.1. Class O-hl and Flat-spectrum YSOs 

The standard method for determining the time spent 
in each YSO Class is to calculate the ratio of the num¬ 
ber of sources in that Class to the number in a reference 
Class, and then multiply by the du r ation of the refer¬ 
ence Class le.g.. IWhking et al.lll9^ lEvans et al.l 120091 : 
iDunham et al.ll20I4D . Implicit in this method are the as¬ 
sumptions that star formation is continuous over at least 
the duration of the reference class and that time is the 
only variable. Regarding the first assumption, we strive 
to average out cloud-to-cloud variations in the recent star 
formation history by averaging over many clouds. Re¬ 
garding the second, our results are best thought of as an 
average over all other relevant variables (e.g., environ¬ 
ment, stellar mass; see lEvans et al.l [20090 . In this sec¬ 
tion, we investigate the effects of different assumptions 
for the reference class and its duration on the timescales 
of Class O-l-I and Flat-spectrum YSOs. _ 

We first follow the method of lEvans et al.l (|2OO90 and 
iDunham et al.l (|20I4D . who adopt Class II as the refer¬ 
ence Class and assume a Class II duration of 2 Myr. 
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[ 5 . 8 ]' - [ 8 . 0 ]' 


Figure 8. Extinction corrected color-color diagrams for the fo ur IRAC bands, wi th the colors and symbols the same as in Figure [4| The 
box in (a) indicates the region identified with Class II sources bv I Allen et al.l l|20Q4ll . whereas in (b) it indicates our proposed modifications 
to this region (see text for details). 



Figure 9. Extinction corrected color-color diagram using three of the four IRAC bands and the first MIPS band, with the colors and 
symbols the same as in Figure!^ The boxes in (a) indicate the regions identified with Class Ill/stellar, Class II, and Class O-hl sources by 
IMuzerolle et all (|20Q4lV going from lower left to upper right. The boxes in (b) indicate our proposed modifications to these regions, now 
with four regions (Class Ill/stellar, Class II, Flat-spectrum, and Class O-l-I, again from lower left to upper right). 
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This duration is taken from analyses of the fraction of 
cluster members with disks versus cluster age for various 
young clusters , which generally show a disk disperal time 
of a few Myr (jWvattl 120081 iRiba s et al.ll 2ni4 1201^ . As 
discussed extensively bv lEvans et al.l (j2009li . this time is 
best considered as a disk half-life rather than an absolute 
duration. The first row of Table 0 presents durations for 
the Class O-l-I and Flat-spectrum Classes calculated fol¬ 
lowing this method, using the extinction corrected num¬ 
bers: 0.52 and 0.34 Myr for Class O-fl and Flat-spectrum 
objects, respectively. 

Although the ages of young stars are highly uncertain, 
recent studies suggest a revised half-life for disk survival 
of 2.5 — 4 Myr on both observational (e.g.. | Kraus et al.l 
l2012t iBell et ani2013t iSoderblom et al.ll2014ll and theo¬ 
retical le.g.. [Alexander fc Armitagell2009D gr ounds. Very 
recent results bv IRibas et al.l ()2014[ 1 and IRibas et'm 
(HnH), using spectroscopically confirmed members of 
various young clusters, show that this duration depends 
on the wavelength where the infrared excess indicative 
of the presence of a disk first appears, ranging from 1.9 
Myr for A < 5 /rm to 3.0 Myr for A < 24 p,m (note that 
these times have been converted from the 1/e times pre¬ 
sented by Ribas et al. to half-life times). As our Spitzer 
data are sensitive to infrared excesses out to 24 ^m, and 
our Spitzer c2d-|-GB data have comparable mid-infrared 
sensitivities to the Ribas et al. Spitzer and WISE data 
used to determine disk lifetimes, their derived lifetimes 
suggest that 3 Myr is a more appropriate choice for the 
duration of the reference Class in this work. The sec¬ 
ond row of Table [7] lists the revised durations: 0.78 and 
0.50 Myr for Class O-l-I and Flat-spectrum objects, re¬ 
spectively. 

While it has become somewhat standard to adopt Class 
II as the reference Class, this choice may result in a logi¬ 
cal inconsistency. The duration of either 2 or 3 Myr based 
on disk fraction vs. cluster age is derived by determining 
the ratio of cluster members with infrared excesses to all 
cluster members (in other words, the ratio of disk sources 
to disk -I- no disk sources). By adopting Class II as the 
reference Class, we omit all of the Class III YSOs that 
still have disks, even though these objects are included 
when other studies (such as those by Ribas et al. ) derive 
disk lifetimes. Since infrared colors indicative of dust are 
the fundamental ingredient in our identifcation of YSOs 
(see 92.2|) . by definition the Class III YSOs identified in 
this work have disks (if they are all in fact YSOs; the issue 
of Class III contamination will be addresse d in the fol¬ 
lowing par agraph). Since the r ecent work bv IRibas et al.l 
(|2014h and IRibas et’aD (I2015D is based on mid-infrared 
data from Spitzer and WISE with comparable sensitivity 
to mid-infrared excesses as our c2d-|-GB data, and since 
we use disk lifetimes derived by Ribas et al. to set the 
duration of the reference Glass, and further given that all 
Glass II and Class III sources identified here have disks, 
we thus argue that it is more appropriate to adopt the 
sum of Class II and Class III as the reference Class than 
Class II alone. Thus, the third row of Table [7] lists the 
revised durations assuming a duration of 3 Myr for the 
sum of Class II and III sources: 0.40 and 0.26 Myr for 
Class 0-fI and Flat-spectrum objects, respectively. 

The one major complication with the above analysis is 
the fact that our Class III sample suffers from 25% - 90% 
contamination, as discussed in 92.31 On the other hand, 


conta mination is not a problem in the IRibas et al.l (|20l4 
1201511 samples since they only consider spectroscopically 
confirmed cluster members. As a result, the above calcu¬ 
lations of Class durations likely overcount the number of 
objects in the reference Class and thus underestimate the 
Class O-fl and Flat-spectrum durations. To account for 
this range of contamination rates, the last two columns 
of Table [7] list the calculated durations of the Class O-l-I 
and Flat-spectrum Classes with the Class II -|- Class III 
reference Class corrected for contamination and a refer¬ 
ence Class duration of 3 Myr. We consider these values 
to be our most accurate and logically consistent lifetimes: 
0.46 — 0.72 Myr for Class 0-(-I sources, and 0.30 — 0.46 
Myr for Flat-spectrum sources. 

4.2. Class 0 YSOs 

In 93.21 we reported that 33.5% of the objects asso¬ 
ciated with dense cores (and thus considered protostars 
according to the definition adopted in this paper) are 
classified as Class 0 by with the remainder classi¬ 
fied as Class I. Under the same assumptions as above 
(continuous star formation and time is the only variable) 
and our definition of a protostar as a YSO associated 
with a dense core, our results imply that 33.5% of the 
total protostellar duration is spent in the Class 0 phase, 
with the remaining 66.5% of the time spent in the Class 
I phase. If we assume that the lifetime of the SED Class 
O-fl objects is a good proxy for the lifetime of the proto¬ 
stellar stage of evolution (see 94.31 below for justification 
of this assumption), the range of SED Class O-fl dura¬ 
tions listed in Table[7] (0.40 — 0.78 Myr) imply resulting 
durations of 0.13 — 0.26 Myr for Class 0 protostars and 
0.27 — 0.52 Myr for Class I protostars (see the last two 
columns of Table [3). If we only consider the last two 
rows of Table [3 (including Class III YSOs in the refer¬ 
ence Class, correcting for estimated contamination rates 
in the Class III sample, and adopting a reference Class 
duration of 3 Myr), the resulting durations are 0.15—0.24 
Myr for Class 0 and 0.31 — 0.48 Myr for Class I. 

As noted in 93.21 we generally find the same ratios of 
SE P Class 0 to SEP Class O-l-I YSOs as in previous stud¬ 
ies (EWans et al.l[20091 : lEnoch et al.ll2009l : iDunham et al.l 
1201411 . While those studies derived a Class 0 duration of 
0.15 — 0.17 Myr, our expanded range of possible dura¬ 
tions is due to our expanded analysis of the contents and 
duration of the reference Class. 

4.3. Lifetimes of Protostars 

As noted in 93.21 the standard picture of low-mass 
star formation resulted from the merger of an empiri¬ 
cal classification system based on the infrared spectral 
indices of YSOs (|Lada fc Wilkind 1 19841 1 with theoretical 
studies o f the stages of the collapse o f a dense, rotat¬ 
ing core (jShdll977hlTerebev et al.lll984[ l. In this picture, 
SED Class O-l-I YSOs are protostars still embedded in 
and accreting from dense cores, SED Class II YSOs are 
T Tauri stars with optically thick, accreting disks, and 
SED Class III YSOs are pre-main sequence stars with lit¬ 
tle or no disks left (all 5pzteer-identified Class III YSOs, 
by definition, still have infrared excesses, thus the pop¬ 
ulation of Class III YSOs without disks is not identified 
in the work presented here). Plat-spectrum YSOs have 
no clear link to a distinct physical stage, and in real¬ 
ity the correlation between observational SED Class and 
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physical Stage is weakened by the e ffects of geometry, ex¬ 
tinction, and accretion history (e.g..lWhitney et al.ll200l: 
Young fc Evansl[200^lRobitaille et al.ll2006l : lCraDsi et al l 

2008HDunham et alJl201fl[l' 

While our study provides the statistics necessary to 
measure durations of the observed SED Classes, it does 
not provide the necessary information to measure dura¬ 
tions of physical Stages. In a sep arate, complementary 
study, iHeiderman fc Evaiii (1201511 surveyed nearly all of 
the SED Class O-l-I and Flat-spectrum YSOs in the c2d 
and GB surveys for emission from the dense gas tracer 
HCO+ (3-2). Using criteria based on the detection and 
strength of this line, they sorted objects into Stage O-l-I 
(protostar) and Stage 11 (disk) sources. They found that 
the Flat-spectrum YSOs are approximately evenly di¬ 
vided between protostar and disk sources, with no evi¬ 
dence that they occupy a physically distinct evolutionary 
stage. They also found that the number of Flat-spectrum 
YSOs classified as protostars is balanced by the number 
of SED Class O-l-I YSOs classified as disks, coinciden¬ 
tally resulting in a nearly identical number of protostars 
(defined in their work based on the presence of dense 
gas traced by HCO+ line emission) and SED Class 0-(-I 
YSOs (326 Class O-l-I YSOs identified here, 335 proto¬ 
stars identified by Heiderman & Evans). Thus, based on 
their results, the range of possible protostellar stage du¬ 
rations is the same as the range of possible SED Class 
O-l-I durations listed in Table [T] 

5. SUMMARY 

In this paper we have presented the full catalog of 
young stellar objects identified in the molecular clouds 
surveyed by the Spitzer Space Telescope “cores to disks” 
(c2d) and “Gould Belt” (GB) Legacy surveys. In the 
case of the c2d c louds, our c a talog updates a previous 
one published bv lEvans et al.l (1200911 . and in the case of 
the GB clouds, our catalog represents the first catalog 
of 6piteer-identified YSOs in the full ensemble of these 
clouds. We summarize our main results as follows: 

• We have used standard techniques developed by 
the c2d project to identify 3239 candidate YSOs 
in the 18 c2d and GB clouds, 2966 of which sur¬ 
vive visual inspection and form our final catalog of 
young stellar objects in the Gould Belt. 

• We compile SEDs for all 2966 YSOs and correct 
these SEDs for foreground extinction. 

• We calculate and tabulate the infrared spectral in¬ 
dex, bolometric luminosity, and bolometric tem¬ 
perature for each object, both with and without ex¬ 
tinction corrections. We also tabulate information 
on whether or not each YSO is considered highly 
likely to be a background AGB star, and whether 
or not each YSO is associated with a dense core as 
traced by submillimeter or millimeter continuum 
emission. 

• After correcting for extinction, we find that 326 out 
of 2966 (11%) are classified as Glass 0-fI, 210 out of 
2966 (7%) are classified as Flat-spectrum, 1248 out 
of 29(36 (42%) are classified as Class II, and 1182 
out of 2966 (40%) are classified as Class III. The 


Class III sample suffers from an overall contamina¬ 
tion rate by background AGB stars between 25% 
and 90%. 

• Adopting standard assumptions, we derive dura¬ 
tions of 0.40 — 0.78 Myr for Glass 0-(-I YSOs and 
0.26 — 0.50 Myr for Flat-spectrum YSOs, with most 
likely durations of 0.46 — 0.72 Myr (Class 0-(-I) and 
0.30 — 0.46 Myr (Flat-spectrum). The ranges en¬ 
compass uncertainties in the adopted assumptions. 

• Including information from submillimeter and mil¬ 
limeter wavelengths, one-third of the Class O-l-I 
sample is classified as Class 0. With this result, 
we calculate durations of 0.13 — 0.26 Myr (Class 0) 
and 0.27 — 0.52 Myr (Class I), with most likely du¬ 
rations of 0.15 — 0.24 Myr (Class 0) and 0.31 — 0.48 
Myr (Class I). As above, the ranges encompass un¬ 
certainties in the adopted assumptions. 

• We revisit infrared color-color diagrams used in the 
literature to classify YSOs and use our large, ex¬ 
tinction corrected sample to propose revisions to 
classification boundaries in these diagrams. 

• The bolometric temperature of a YSO surrounded 
by a circumstellar disk but no longer embedded in 
a dense core is primarily determined by the spec¬ 
tral type of the YSO and the availability of optical 
photometry, and not by the properties of the cir¬ 
cumstellar material. As such, the bolometric tem¬ 
perature should not be used to distinguish between 
Class II and Class III YSOs. 

The YSO catalog presented here represents the most 
complete catalog of star formation within 500 pc of the 
Sun yet assembled. It includes nearly every major site 
of star formation in the solar neighborhood except for 
the Taurus and Orion molecular clouds, each of which 
was the subject of its own, dedicated Spitzer Legacy sur¬ 
vey. While this work represents a major step forward 
in assembling a complete and reliable catalog of Young 
Stellar Objects in the Gould Belt, future work to further 
refine the YSO catalog presented here is needed on the 
following three fronts: (1) combining the results of the 
c2d-|-GB, Taurus, and Orion surveys using uniform data 
reduction and analysis techniques, (2) identifying and re¬ 
moving background AGB stars currently contaminating 
the sample, and (3) revisiting the associations between 
Spitzer sources and dense cores traced by submillimeter 
and millimeter detections once the full results from the 
Herschel and SCUBA-2 Gould Belt surveys are available. 
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APPENDIX 

A. SPITZER ASTRONOMICAL OBSERVATION REQUESTS (AORS) FOR UNPUBLISHED CLOUDS 

Three clouds listed in Table [1] do not have published data references (Chamaeleon I, Cha maeleon HI. and Musca) . 
Additionally, only a small piece of Aquila was published in the reference listed for that cloud (|Gutermuth et al]l2008fl . 
Consequently, for these four clouds there are no published lists of the Spitzer Astronomical Observation Requests 
(AORs) that targeted them. We thus list their AORs in Table El 

B. ROUNDING ARTIFACTS IN COLOR-COLOR DIAGRAMS 

Inspection of Figure [2] reveals a vertical line of sources at [3.6]-[4.5] ^ 0.5. If an object has identical flux densities 
at two Spitzer wavelengths, the magnitude difference between these two wavelengths is then simply equal to: 

mi - TO 2 = 2.5 log , (Bl) 

jjzmi 

where ffzmi and f fzm^ are the fluxes for zero magnitude at each wavelength. For 3.6 and 4.5 /im, this gives 
[3.6]-[4.5] = 0.473. 

There are 147 objects plotted in Figure E] with identical flux densities at 3.6 and 4.5 ^m, resulting in 147 objects 
plotted at the same value on the x-axis and thus giving rise to the vertical line of sources. We note that this is 
an artifact introducing by rounding all flux densities to two significant digits. To demonstrate this point, Figure (TU] 
plots the same color-color diagram as Figure El except now using the original data before rounding. As evident from 
inspection of this figure, the vertical line of sources at [3.6]-[4.5] ^ 0.5 has now disappeared. 

The same artifact is also present in Figure El in the [3.6]-[24] color. For 3.6 and 24 /im, identical flux densities 
result in a magnitude difference of [3.6]-[24] = 3.959. There are only 23 objects plotted in Figure El with identical 
flux densities at 3.6 and 24 /im, so while the vertical line of sources at [3.6]-[24] ^ 4 is present, it is not as apparent. 
Similar artifacts also affect the colors plotted in Figures E] and El using the extinction corrected photometry. For the 
four colors plotted in these two figures, between 38 and 95 objects with identical flux densities are plotted. While 
corresponding vertical and horizontal lines of sources are present, none are apparent because they all occur at colors 
where the source densities are too high for these small subsets of sources to be visible. These rounding artifacts have 
no effect on the results presented in this paper. 
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[ 3 . 6 ]-[ 4 . 5 ] 


Figure 10. Same as Figure[2] ([3.6]-[24] versus [3.6]-[4.5] for all 2966 objects in our final YSO sample, using the observed photometry), 
except plotted using the original data before rounding. The vertical line of sources at [3.6]-[4.5] ~ 0.5 has now disappeared as YSOs with 
similar flux densities have not now all been rounded to identical flux densities. 


C. BOLOMETRIC TEMPERATURES OF CLASS II AND III SOURCES 


As noted in 93.11 the bolometric temperatures of our Class 11 and 111 YSOs are generally tightly clustered into 
two narrow ranges (see Figures 0] and [7]). To understand the origins of these effects, we construct simple con¬ 
tinuum radiative transfer models of circumstellar disks using the three-dimensional Monte Carlo radiative transfer 
package RADMC-3D^®. We assume an axisymmetric physical structure and use RADMC-3D to first calculate the two- 
dimensional dust temperature profile through the disks and then generate full model SEDs from optical to millimeter 
wavelengths. 

We adopt a simple, parameterized disk density structure following a power-law in the radial coordinate and a 
Gaussian in the vertical coordinate. The exact density profile is given as follows: 


Pdiski^: ■2') — Po 




(Cl) 


where z is the distance above the midplane (z = rcosO, with r and 9 the usual radial and zenith angle spherical 
coordinates), s is the distance in the midplane from the origin (s = \/t^ — z^), and po is the density at the reference 
midplane distance sq, normalized such that the total disk mass is correct (see below). The quantity Hg is the disk 

scale height and is given by Hs = Ho (j") ^ where Ho is the scale height at So- The parameter /3, which describes 

how the scale height changes with s and sets the flaring of the disk, is treated as a free parameter (see below). We set 
Ho = 10 AU at So = 100 AU. The disk surface density profile, calculated by integrating Eg nation I C1 1 over the vertical 
coordinate z, has a radial power-law index of E(s) oc where p = a — (3. For these models, we hold a fixed at 2.25 
and the disk outer radius fixed at 100 AU. 

To consider a diverse set of disks, we treat the disk inner radius, flaring parameter /3, and total disk mass as free 
parameters. Allowed values for each parameter are as follows: 0.1, 1, 10, 20, and 50 AU for the disk inner radius; 1.0, 
1.25, 1.5, 1.75 for /3, and 0.0005, 0.001, 0.005, 0.01, and 0.05 Mq for the disk mass. We thus consider 100 total models 
and generate SEDs at 9 inclinations^^, starting at z = 5° and increa sing in steps of 10°, resulting in 900 total model 
SEDs. We adopt the dust opacities of lOssenkopf fc Hennind (1199411 appropriate for thin ice mantles after 10® yr of 
coagulation at a gas density of 10® cm“® (OHS dust). The radiative transfer is performed with no external heating 
by the interstellar radiation field, but isotropic scattering off dust grains is included. The input stellar spectrum is 
assumed to be a simple blackbody with L = 1 Lq and Teff of eith er 4050 K or 3050 K (corr esponding approximately 
to K7 and M5 spectral type pre-main sequence stars, respectively: iPecaut fc Mamaiekll2M^ . Finally, for each of the 
1800 total model SEDs (two stellar models each surrounded by 100 different circumstellar disk models, each viewed at 


Available at: http://www.ita.uni-heidelberg.de/~dullemond/software/radmc-3d/ 

An inclination of i = 0° corresponds to a pole-on (face-on) system, while an inclination of i = 90° corresponds to an edge-on system. 
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Bolometric Temperature (K) 


Figure 11. Histograms showing the distribution of for the models described in the text. Left: The distributions of for the K7 
disk models, calculated both with (shaded) and without (dashed) optical photometry included in the SEDs. Right: The distributions of 
'^bol (shaded) and M5 (dashed) disk models, in both cases calculated with optical photometry included in the SEDs. 

9 inclinations), we calculate two values of the bolometric temperature: one with the full SED from optical to millimeter 
wavelengths, and another excluding all optical wavelengths (all wavelengths shortward of 1.25 /im). 

The left panel of Figure [TT] shows the distributions of for the K7 disk models, calculated both with and without 
optical photometry included in the SEDs. Since the radiative transfer model includes no foreground extinction, these 
values are equivalent to our observed values after correcting for extinction. When optical photometry is included in 
the SEDs, the distribution of exhibits a narrow peak around 3500 K with a width of only a few hundred K, plus a 
long tail to lower values corresponding to SEDs viewed at edge-on inclinations. A similar behavior is exhibited when 
optical photometry is left out of the SEDs, except both the peak and long tail are shifted to lower values of This 
behavior is expected since the inclusion of shorter wavelength data shifts the flux-weighted mean frequency to higher 
values. 

The right panel of Figure [TT] shows the distributions of for the K7 and M5 disk models, in both cases calculated 
with optical photometry included in the SEDs. Each set of models exhibits narrow peaks, with the K7 model peak 
located at higher values than the M5 model. Even for the large range of disk properties considered here, it is evident 
from this Figure that the measured for YSOs with disks is mostly determined by the spectral type of the central 
source. The disk properties introduce only a very narrow spread except for models viewed at nearly edge-on inclinations. 

The tight clustering into two narrow ranges of for most of our Class II and III YSOs, as evident from Figures ID 
and]?] is naturally explained by these results. A K7 spectral type is assumed for most sources when deriving extinction 
corrections, explaining the tight clustering, and optical photometry is only available for a subset of our YSOs (see 
>12.41) . explaining the bifurcation into two separate ranges of One consequence of t his situation is th at, even 

though a boundary in bolometric temperature between Class II and Class III is defined bv iChen et al.l (|I995l) . is 
not a good discriminator between these two Classes. 

D. TABLES OF ANCILLARY PHOTOMETRY COLLECTED FROM THE LITERATURE 

Here we tabulate all ancillary photometry used to construct YSO SEDs, where we define ancillary to mean all 
photometry beyond the 2MASS and Spitzer photometry measured, band-merged, and cataloged by the c2d pipeline. 
This definition of ancillary includes Spitzer 160 pm photometry. All flux densities and flux density uncertainties are 
listed in mJy and rounded to two significant digits. 

Tables [9] and [T0| list, for each YSO, the same running index as in Table |2l followed by flux density and flux density 
uncertainty pairs for all photometry compiled between 0.36 — 0.96 pm. Table |9| lists observed values whereas Table ITIil 
lists extinction corrected values. 

Tables [TT] and [T2| list, for each YSO, the same running index as in Table [2] followed by flux density and flux density 
uncertainty pairs for all ancillary photometry compiled between 3.4 — 22 pm. Table [TT] lists observed values whereas 
Table [T^ lists extinction corrected values. 

Tables [13] and [14] list, for each YSO, the same running index as in Table [2] followed by flux density and flux density 
uncertainty pairs for all ancillary photometry compiled between 3.4 — 22 pm. Table [T3l lists observed values whereas 
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Table [HI lists extinction corrected values. 


Table 1 

Molecular Clouds Surveyed by the c2d and GB Surveys 


19 


Cloud 

Survey 

Distance 

(pc) 

Distance 

Reference^ 

Area 

(pc^) 

Mass 

(Mq) 

Data 

Reference(s)^ 

Aquila 

GB 

260 

1 

179.0 

24400 ± 3000 

1 

Auriga/CMC 

GB 

450 

2 

52.4 

4844 ± 1152 

2 

Cepheus 

GB 

200-325'= 

3 

38.0 

2610 ± 170 

3 

Chamaeleon I 

GB 

150 

4 

9.4 

857 ± 210 


Chamaeleon II 

c2d 

178 

5 

9.2 

637 ± 300 

4, 5, 6 

Chamaeleon III 

GB 

150 

4 

28.0 

1330 ± 390 


Corona Australis 

GB 

130 

6 

3.0 

279 ±110 

7 

IC5146 

GB 

950 

7 

149.0 

8550 ± 2170 

8 

Lupus I 

c2d 

150 

8 

8.9 

513 ± 310 

9, 10 

Lupus III 

c2d 

200 

8 

15.4 

912 ± 320 

9, 10 

Lupus IV 

c2d 

150 

8 

2.5 

189 ± 95 

9, 10 

Lupus V 

GB 

150 

8 

11.7 

705 ± 220 

11 

Lupus VI 

GB 

150 

8 

6.7 

455 ± 140 

11 

Musca 

GB 

160 

9 

6.8 

335 ± 110 


Ophiuchus 

c2d 

125 

10 

29.6 

3120 ± 1800 

12 

Ophiuchus North 

GB 

130 

11 

7.3 

621 ± 17 

13 

Perseus 

c2d 

250 

12 

73.2 

6590 ± 3600 

14, 15 

Serpens 

c2d 

429 

13 

17.0 

2340 ± 640 

16, 17, 18 


^ Reference for the distances quoted in this work: (1) iMaury et al.| 1 120111) : (2) ILada et al.| 1 120091) : (3) IKirk et al.| | |2009t) : (4) IBelloche et all 
| |2011bD : ISl IWhittet et al.| | |1997D ; f6) INeuhauser &: Forbrichl | |2008l) : f7) |Harvev et al.| | |200SD : f8) IComer6n| | |200gD : f9) |Knude &: Hod | |199SD : (10) 
IWilking et all | |200SD : fll) lde Geus et al.| | |1989D : (12) |Enoch et all J2006I) : (13) IDzib et all | |2011|) . 

^ References presenting the Spitzer IRAC and MIPS observations: (1) IGutermuth et all | |2008t) : (2) IBroekhoven-Fiene et all (|2014[) : (3) 
IKirk et all | |2009l) : (4) lYoung et all l |2005l) : (5) IPorras et al.| l |2007t) : (6) lAlcala et al.| | |200gD : (7) IPeterson et all | |2011D : (8) IHarvey et all 1 )20081) : 
(9) IChapman et all | |2007l) : (10) IMerin et al.| | |200g|) : (11) ISpezzi et al.| l |201lD : (12) IPadgett et all | |200gD : (13) IHatchell et al.] | |2012D : (14) 
IJdrgensen et all ( |200^ : IRebull et all | |2007t) : (16) IHarvey et al.| | |2006l) : (17): [TTarvey et all l |2007b|) : (18) IHarvey et all | |2007a|) - 
^ Different regions within Cepheus are located at different distances; see Kirk et al. (2009) for details. 


Table 2 

Properties of the c2d+GB YSOs 


Index 

Cloud 

Spitzer 

Source Name 
(SSTc2d or SSTgb +) 

Av 

(mag) 

a 

Observed 

T'bol 

(K) 

Lbol 

(L©) 

Extinction Corrected 
'^bol' Lhoi' 

«' (K) (Lo) 

Likely 

AGB? Core? 

1 

Aquila 

J180144.8-044941 

8.2 

-2.16 

1700 

3.0 

-2.52 

2000 

8.2 

N 

2 

Aquila 

J180145.8-045902 

8.0 

-2.14 

1700 

1.4 

-2.56 

2000 

3.7 

N 

3 

Aquila 

J180154.2-043753 

12.4 

-2.57 

2100 

0.52 

-3.22 

2600 

4.7 

Y 

4 

Aquila 

J180229.2-051552 

8.3 

-2.03 

1700 

3.6 

-2.38 

2000 

9.8 

N 

5 

Aquila 

J180302.0-044233 

12.4 

-2.45 

2100 

0.11 

-3.15 

2500 

0.92 

N 

6 

Aquila 

J180308.0-043523 

8.0 

-2.18 

1700 

0.24 

-2.52 

2000 

0.64 

Y 

7 

Aquila 

J180313.8-042845 

6.9 

-2.22 

1800 

0.72 

-2.53 

2000 

1.7 

Y 

8 

Aquila 

J180423.9-042248 

9.0 

-2.15 

1700 

2.2 

-2.44 

2000 

6.3 

Y 

9 

Aquila 

J180435.7-044700 

8.1 

-2.26 

1700 

1.9 

-2.65 

2000 

5.3 

Y 

10 

Aquila 

J180444.5-043706 

12.4 

0.53 

540 

0.055 

-0.10 

1000 

0.11 

N 


The full table is published online in a machine-readable format. Only the first ten lines are shown here. 
















































































































Table 3 

Observed Spitzer and 2MASS Flux Densities, in mjy, of the c2d+GB YSOs 


to 

o 


Index 

J 

1.25 pva 

AJ 

1.25 

H 

1.65 pm 

AH 

1.65 pm 

Ks 

2.17 /^m 

AKs 
2.17 fim 

IRAC 
3.6 pm 

AIRAC 
3.6 pm 

IRAC 

4.5 

AIRAC 

4.5 

IRAC 
5.8 pm 

AIRAC 
5.8 pm 

IRAC 
8.0 /^m 

AIRAC 

8.0 /^m 

MIPS 
24 fim 

AMIPS 
24 fim 

MIPS 
70 /im 

AMIPS 
70 /^m 

1 

290 

7.3 

820 

33 

1100 

21 

620 

36 

420 

24 

350 

17 

220 

11 

67 

6.2 



2 

140 

3.3 

370 

14 

480 

9.3 

300 

15 

160 

7.7 

150 

6.9 

92 

4.3 

36 

3.3 



3 

220 

5.2 

180 

3.4 

120 

2.3 

53 

2.8 

34 

1.7 

23 

1.1 

13 

0.65 

3.2 

0.38 

36 

9.0 

4 

340 

9.9 

890 

36 

1200 

30 

800 

41 

530 

26 

490 

24 

330 

16 

100 

9.6 



5 

37 

1.0 

40 

1.1 

30 

0.83 

13 

0.64 

8.8 

0.42 

5.9 

0.28 

3.5 

0.17 

1.5 

0.22 



6 

23 

0.62 

67 

1.6 

81 

1.7 

49 

2.5 

33 

1.6 

26 

1.3 

19 

0.88 

4.7 

0.48 



7 

88 

1.9 

200 

4.8 

250 

4.9 

140 

7.4 

96 

4.7 

77 

3.7 

50 

2.4 

14 

1.3 



8 

180 

4.3 

540 

18 

750 

18 

510 

25 

340 

17 

310 

15 

200 

9.5 

48 

4.5 



9 

200 

4.8 

520 

16 

700 

16 

410 

21 

230 

11 

200 

9.6 

120 

5.8 

40 

3.7 



10 

0.62 

0.050 

2.2 

0.090 

3.6 

0.14 

4.7 

0.23 

6.9 

0.33 

12 

0.57 

27 

1.3 

67 

6.2 

110 

16 


a 


The full table is published online in a machine-readable format. Only the first ten lines are shown here. 










Table 4 

Extinction Corrected Spitzer and 2MASS Flux Densities, in mjy, of the 
c2d+GB YSOs 


Index 

J 

1.25 fj.m 

AJ 

1.25 

H 

1.25 /im 

AH 

1.65 fim 

Ks 

2.17 Atm 

AKs 
2.17 /tm 

IRAC 
3.6 fim 

AIRAC 
3.6 fim 

IRAC 

4.5 

AIRAC 

4.5 

IRAC 
5.8 fim 

AIRAC 
5.8 fim 

IRAC 

8.0 Atm 

AIRAC 

8.0 Atm 

MIPS 
24 Atm 

AMIPS 
24 Atm 

MIPS 
70 Atm 

AMIPS 
70 Atm 

1 

2200 

55 

2800 

120 

2400 

49 

1000 

59 

630 

35 

490 

23 

320 

16 

82 

7.6 



2 

1000 

24 

1300 

46 

1100 

21 

490 

25 

240 

12 

200 

9.6 

130 

6.3 

43 

4.0 



3 

4700 

no 

1200 

22 

430 

8.3 

no 

5.9 

63 

3.1 

39 

1.8 

24 

1.2 

4.4 

0.51 

38 

9.4 

4 

2600 

78 

3200 

130 

2900 

70 

1300 

68 

800 

40 

680 

34 

490 

23 

130 

12 



5 

810 

22 

260 

7.2 

no 

3.0 

28 

1.4 

16 

0.78 

9.7 

0.46 

6.3 

0.30 

2.0 

0.30 



6 

170 

4.5 

220 

5.4 

190 

3.9 

79 

4.0 

49 

2.3 

37 

1.7 

27 

1.3 

5.7 

0.58 



7 

480 

11 

570 

14 

510 

9.9 

220 

11 

140 

6.6 

100 

4.9 

69 

3.3 

17 

1.6 



8 

1700 

40 

2100 

69 

1900 

45 

880 

44 

530 

27 

450 

21 

300 

14 

60 

5.6 



9 

1500 

35 

1800 

53 

1600 

36 

680 

34 

340 

17 

270 

13 

170 

8.4 

48 

4.5 



10 

13 

1.1 

14 

0.59 

13 

0.50 

10 

0.49 

13 

0.61 

20 

0.95 

48 

2.2 

90 

8.4 

no 

17 


The full table is published online in a machine-readable format. Only the first ten lines are shown here. 
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Table 5 

Number of YSOs by Cloud and Class 


Extinction Corrected, All Sources Extinction Corrected, Excluding Likely AGB Stars 


Cloud 

O+I 

Flat 

II 

III 

Total 

O+I 

Flat 

II 

III 

Total 

Aquila 

83 

65 

330 

841 

1319 

83 

65 

327 

634 

1109 

Auriga/CMC 

35 

8 

73 

17 

133 

35 

8 

73 

15 

131 

Cepheus 

18 

11 

61 

13 

103 

18 

11 

61 

12 

102 

Chamaeleon I 

4 

5 

62 

15 

86 

4 

5 

62 

14 

85 

Chamaeleon II 

2 

0 

19 

5 

26 

2 

0 

19 

4 

25 

Chamaeleon III 

1 

0 

0 

3 

4 

1 

0 

0 

0 

1 

Corona Australis 

9 

6 

22 

17 

54 

9 

6 

22 

16 

53 

IC5146 

28 

10 

79 

15 

132 

28 

10 

79 

14 

131 

Lupus I 

1 

2 

8 

2 

13 

1 

2 

8 

0 

11 

Lupus III 

2 

5 

38 

24 

69 

2 

5 

38 

14 

59 

Lupus IV 

1 

1 

3 

7 

12 

1 

1 

3 

5 

10 

Lupus V 

0 

0 

7 

36 

43 

0 

0 

7 

25 

32 

Lupus VI 

1 

0 

2 

42 

45 

1 

0 

2 

25 

28 

Musca 

1 

0 

1 

11 

13 

1 

0 

1 

4 

6 

Ophiuchus 

28 

43 

174 

47 

292 

28 

43 

174 

32 

277 

Ophiuchus North 

2 

1 

3 

4 

10 

2 

1 

3 

3 

9 

Perseus 

76 

35 

235 

39 

385 

76 

35 

235 

31 

377 

Serpens 

34 

18 

131 

44 

227 

34 

18 

131 

39 

222 

Total 

326 

210 

1248 

1182 

2966 

326 

210 

1245 

887 

2668 


Table 6 

Classification Via [3.6]^ — [5.8]^ versus [8.0]^ — [24]' 


Quantity 

Class Ill/stellar 

Class II 

Flat 

Class O+I 

Original x^in 

- 0.5 

2.0 


3.50 

Original Xmax 

1.5 

5.0 



Original ymin 

- 0.5 

0.3 


1.50 

Original ymax 

0.5 

1.5 



Original Pi^ 

65.2% 

77.9% 


66.0% 

Original P 2 ^ 

99.5% 

84.2% 


84.1% 

Revised Xrr.ir. 

-0.5 

1.5 

2.5 

2.5 

Revised Xmax 

4.0 

5.0 

4.5 


Revised y^in 

-0.5 

0.4 

1.1 

1.9 

Revised ymax 

0.4 

1.1 

1.8 


Revised Pi^ 

83.2% 

79.2% 

63.2% 

65.6% 

Revised P 2 *’ 

92.2% 

91.8% 

59.1% 

90.0% 


The probability that a YSO of the desired Class is located within the region. 
The probability that a YSO located within the region has the desired Class. 


Table 7 

Timescales for Young Stellar Objects 


Reference Class 

(Myr) 

(Myr) 

(Myr) 

(Myr) 

(Myr) 

All Class II 

2.0 

0.52 

0.34 

0.17 

0.35 

All Class II 

3.0 

0.78 

0.50 

0.26 

0.52 

All Class II + All Class III 

3.0 

0.40 

0.26 

0.13 

0.27 

All Class II + 75% of Class III 

3.0 

0.46 

0.30 

0.15 

0.31 

All Class II + 10% of Class III 

3.0 

0.72 

0.46 

0.24 

0.48 


^ Assumed duration of the reference Class. 

^ Calculated duration of Class O+I sources (Tref-jvr^)- 

^'ref 

Calculated duration of Flat-spectrum sources (rref jv )• 

^ Calculated duration of Class 0 sources (0.335 X tq+i); see ^4.2l for details. 
^ Calculated duration of Class I sources (0.665 X tq+i); see )j4.2| for details. 


Table 8 

Spitzer AORs for Unpublished Clouds 


Cloud 

MIPS Target Name 

MIPS AOR 

IRAC Target Name 

IRAC AOR 

Aquila 

SERAQUJl.l 

SERAQUJ1.2 

19968768 

19968256 

20001536 

SERAQUJl.l 

20008448 

20008704 

20008960 
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Table 8 — Continued 


Cloud 

MIPS Target Name 

MIPS AOR 

IRAC Target Name 

IRAC AOR 



20001024 


19970816 


seraquj:.3 

19970560 


19971072 



19970304 


19971328 


seraquj:.4 

20002048 

seraquj:.2 

20005120 



20001792 

seraquj:.3 

20005376 

20005632 

19964672 

19965440 

19965952 

20002304 

20002560 

20003072 

19957248 

19958016 

19958272 




seraquj:.4 

19997952 

19998464 

19998720 

19999488 

20018176 

20018688 

20018944 

20019200 




seraquj:.5 

14511104 




SERAQUJ:.6 

27042304 

27042560 


SERAQUJIE 

14510336 

SERAQUJIE 

19991552 



19966464 


19988480 



19966976 


19991296 

20009728 

19987968 

20009216 


SERAQU.SW.1 

24195328 

SERAQU.SW 

19972352 



24195584 


19971840 


SERAQU.SW.2 

24195840 


19972608 



24196096 


19986944 

19986688 

19987200 


SERAQU.W.l 

19974656 

SERAQU.W.la 

19970048 



19974400 

SERAQU.W.lb 

19985664 

19962368 

19980800 


SERAQU.W.2 

20006912 

SERAQU.W.2 

19957760 



20006656 


19956992 

19975424 

19974912 

Chamaeleon I 

CHAJ 

03661312 

CHAJ.1 

19986432 



03962112 


20006400 



19979264 

CHA.I.2 

03960320 



20011008 

CHA.I.3 

20015104 



19978496 


20014592 



20010240 

CHA.I.4 

03651328 




CHA.I.5 

19992832 

20012800 

Chamaeleon III 

CHAJII.l 

19980032 

CHAJII.l 

19995392 



19979776 


20015872 


CHAJII.2 

20012544 

CHAJII.2a 

19989248 



19977472 


19989760 



20012032 


20007168 



19976448 

CHAJII.2b 

20007424 

19985152 

20003328 




CHAJII.2C 

19961344 

19977216 

19984896 

19961600 

19998208 

20002816 

27040256 

27040512 

27041536 

27042048 

Musca 

MUSCA 

20004096 

MUSCA.1 

27040768 



19973120 


27041024 



20003840 

MUSCA.2 

27043072 



19972864 

MUSCA.3 

27043328 

20009472 
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Table 8 — Continued 


Cloud MIPS Target Name MIPS AOR IRAC Target Name IRAC AOR 


19975680 

19976704 

20009984 

19975936 

19977728 



CDOO-qOiCnji^COtO 


Table 9 

Observed 0.36 — 0.96 /rm Flux Densities, in mjy, of the c2d+GB YSOs 


u 

AU 

B 

AB 

V 

AV 

R 

AR 

Hal2 

AHal2 

I 

AI 

m914 

Am914 z Az 

Index 0.36 /im 

0.36 fj.m 

0.44 /tm 

0.44 fim 

0.55 

0.55 fim 

0.64 

0.64 

0.665 /rm 

0.665 /tm 

0.79 /rm 

0.79 /tm 

0.915 /rm 

0.915 /rm 0.96 /rm 0.96 /tm 


1 


10 

^ The full table is published online in a machine-readable format. Only the first ten lines are shown here. 
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cDOo-qaicnji^coto 


Table 10 

Extinction Corrected 0.36 — 0.96 /im Flux Densities, in mjy, of the c2d+GB cti 

YSOs 


u 

AU 

B 

AB 

V 

Ay 

R 

AR 

Hal2 

AHal2 

I 

AI 

m914 

Am914 z Az 

Index 0.36 /im 

0.36 /rm 

0.44 fim 

0.44 fim 

0.55 

0.55 fim 

0.64 

0.64 

0.665 fim 

0.665 fim 

0.79 fim 

0.79 fim 

0.915 fim 

0.915 fim 0.96 fim 0.96 fim 


1 


10 

^ The full table is published online in a machine-readable format. Only the first ten lines are shown here. 























































Table 11 

Observed 3.4 — 22 /rm Flux Densities, in mjy, of the c2d+GB YSOs 


27 


Index 

L AL M AM 

3.4 fj.m 3.4 fj,m 5.0 /im 5.0 /i.m 

ISO 
6.7 /im 

Also 
6.7 /im 

WISE 

12 fj,m 

AWISE 

12 fim 

ISO 

14 fj,m 

Also 

14 fim 

WISE 

22 fim 

AWISE 

22 fim 

1 




160 

2.3 



87 

2.3 

2 




70 

1.0 



47 

1.6 

3 




7.6 

0.21 



3.7 

0.99 

4 




240 

3.3 



120 

2.9 

5 




1.7 

0.16 





6 




14 

0.28 



8.2 

0.85 

7 




33 

0.51 



19 

1.1 

8 




120 

1.5 



63 

2.1 

9 




81 

1.2 



50 

1.7 

10 




27 

0.47 



65 

2.2 

The full table is published online in a machine-readable format. Only the first ten lines are shown here. 



Extinction Corrected 3.4 — 22 

Table 12 

lira Flux Densities, in mJy, of the c2d+GB 
YSOs 



Index 

L AL M AM ISO 

3.4 /im 3.4 /i.m 5.0 /im 5.0 /i.m 6.7 /i.m 

Also 

6.7 /im 

WISE AWISE 

12 fj,m 12 /im 

ISO Also 
14 fj,m 14 fim 

WISE AWISE 

22 fj,m 22 /im 

1 



230 

3.4 


110 

2.9 

2 



100 

1.5 


59 

2.0 

3 



14 

0.39 


5.2 

1.4 

4 



360 

4.9 


160 

3.6 

5 



3.1 

0.29 




6 



21 

0.42 


10 

1.1 

7 



46 

0.72 


23 

1.3 

8 



190 

2.4 


81 

2.7 

9 



120 

1.8 


62 

2.1 

10 



49 

0.87 


92 

3.2 


The full table is published online in a machine-readable format. Only the first ten lines are shown here. 




















































CDOO-qOiCnji^Coto 


Table 13 

Observed 160 — 1300 /rm Flux Densities, in mjy, of the c2d+GB YSOs 


to 

00 


MIPS AMIPS 

Index 160 /rm 160 /rm 350 /rm A350 /rm 450 /rm A450 /rm 850 /rm A850 /rm 1100 /rm AllOO /rm 1200 /rm A1200 /rm 1300 /rm A1300 /rm 


10 


The full table is published online in a machine-readable format. Only the first ten lines are shown here. 













































Table 14 

Extinction Corrected 160 — 1300 fim Flux Densities, in mjy, of the c2d+GB 

YSOs 


MIPS AMIPS 

Index 160 fim 160 fim 350 fim A350 fim 450 /rm A450 fim 850 fim A850 fim 1100 /rm AllOO fim 1200 fim A1200 /rm 1300 fim A1300 fim 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 


The full table is published online in a machine-readable format. Only the first ten lines are shown here. 













































